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ABSTRACT 


A methodology has been developed to predict the failure of stiffened 
composite panels by stiffener- skin separation. The methodology is applicable to 
curved or flat panels under combined uniaxial compression and in-plane shear loads. 
The analytical crux of the methodology are two analysis packages: (1) PACL1 which 

predicts the stress and displacement fields in a stiffened postbuckled panel and (2) 
WEBSTER, which predicts the local interfacial (shear and normal) stresses at the skin 
and stiffener flange junction. PACLl is a Rayleigh-Ritz analysis of curved composite 
panels loaded into the postbuckling range. WEBSTER is a two-dimensional elasticity 
analysis for the skin- stiffener interface. This paper presents results of the pre- 
dictive methodology. The entire analysis can be performed on a personal computer. 


INTRODUCTION 


Recent studies have indicated that the most common mode of failure for 
stiffened composite panels loaded into the postbuckling range is stiffener/skin 
separation. An assessment of the current postbuckled stiffened panel design, analy- 
sis, and application technology (References 1 and 2) shows that, for both static and 
fatigue loading, stiffener/skin separation consistently occurs at loads below other 
competing failure modes, such as skin compression failure and delamination. Because 
stif fener/skin separation is not accounted for in the existing design methodology, 
the application of existing postbuckling methodology to the design of advanced com- 
posite panels has resulted in unconservative designs. In order to fully realize the 
weight saving potential of postbuckled designs, this failure mode must be accounted 
for . 


^ This work was performed under NASA/Northrop Contract NASI- 18842, entitled "Innova- 
tive Composite Fuselage Structures." 
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Several attempts have been made in the last decade to predict the 
stiffener/skin separation of postbuckled, stiffened composite panels. These include 
semi-empirical models of References 1 and 3, strength of material models in Referenc- 
es 4, 5, and 6, and elasticity models of References 7 and 8. In Reference 1, an 
empirical equation is derived by analogy with the crippling data for plates with one 
edge simply supported and one edge free. It is hypothesized that when the panel web 
strain reaches the crippling strain, the interfacial stresses become high enough to 
cause failure. The model proposed in Reference 3 is based on the hypothesis that 
skin/stiffener separation is caused by the interface normal stress Induced by 
torsional stress in the stiffener. The average torsional stress in the stiffener, 
obtained from postbuckling analysis method of Reference 9 is empirically related to 
the bondline normal stress. 

In Reference 4, a stiffener/skin delamination model was developed to 
determine the out-of-plane stresses in postbuckled stiffened shear panels. The 
physical behavior is simulated by a beam aligned parallel to the diagonal buckles 
through the point of maximum deflection. The Rayleigh-Ritz method was used to obtain 
the interfacial stresses. A similar model based on the large deflection beam theory 
was developed in References 5 and 6. 

A three-dimensional analysis, based on the principle of virtual work was 
developed in Reference 7. The problem was formulated based on plate theory; continu- 
ity was enforced through the adhesive interface deformation. The problem was solved 
by the Galerkin method. Three out-of -plane stress components were obtained. A two- 
dimensional local elasticity model was used In Reference 8 to investigate the state 
of stress at the interface between a composite skin and the flanges of a cocured or 
secondarily bonded stiffened panel. 

The preceding investigations examined the local response at the skin/stiff- 
ener interface. No attempts were made to relate the remote loading to the local 
response. In this paper, a methodology is developed that integrates the postbuckling 
analysis method of Reference 10 and the out-of-plane stress analysis model of Refer- 
ences 8 and 11 into a single analysis package. The analytical approach and numeri- 
cal results are presented in the following sections. The entire analysis can be per- 
formed on a personal computer. 


ANALYSIS METHODOLOGY 


The crux of the present methodology are two analysis packages: (1) PACL1 , 

which predicts the stress and displacement fields in a stiffened postbuckled panel 
and (2) WEBSTER, which predicts the local interfacial (shear and normal) stresses at 
the skin and stiffener flange junction. PACL1 (Reference 10) is a Rayleigh-Ritz 
analysis of curved composite panels loaded into the postbuckling range under axial 
compression and shear loads. WEBSTER (Reference 11) is a two-dimensional elastici ty 
analysis for the skin/stiffener interface. These analysis methods and the integrated 
methodology are briefly discussed in the following paragraphs. 

PACL1 . The energy approach was used to formulate the postbuckling problem 
in the PACL1 analysis package. The problem was formulated for a stiffened, curved 
anisotropic laminated plate. The laminate was assumed to be balanced and symmetric. 

A small imperfection in the lateral displacement was also included in the formula- 
tion. The panel geometry and the coordinate system are shown in Figure 1. This 
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figure also shows the relationship between the overall postbuckling structural con- 
figuration and the panel geometry used in the analysis. Since the adjacent bays are 
assumed to deform in an identical fashion, a single bay was analyzed. Figure 1 shows 
that the material properties for the skin are A*fj and D*fj , where A*fj (A*ff, &*12 > 
A* 22 » anc * A* 56 > is the skin stiffness matrix and D*fj (D*ff, D *12 » D *16> D *22> D *26> 
and D*6$) is the skin rigidity matrix. The material axes 1 and 2 are assumed to 
coincide with the panel geometry coordinate axes x and y, respectively. The panel, 
with length a and width b, is bounded by stringers along the straight edges and 
frames or rings along the curved edges. The cross-sectional area, Young's modulus 
and moment of inertia of the stringers are A s , E s , and I s , respectively. Those of 
the frames are Af, Ef, and If. The radius of the curved panel is R. 

The energy expressions are written in terms of the displacement components 
u, v, and w in the x, y, and z directions, respectively. The panel is assumed fixed 
along x - o and subjected to a system of combined compression (N x ) and shear (N X y) 
load along the edge x * a. The out-of-plane displacement (w) is assumed to be zero 
along all four edges. 

The total potential energy, II, is the sum of the strain energy stored in the 
skin, U w , in the stringers, U s , in the frames, Uf, and the potential of the external 
loads, fl, and is written as 

n = U w + U s + U f + 0 (1) 

The strain energy in the skin for an anisotropic plate with A ^ » A 26 = 0 
is given by 

* 2 * * 2 * 2 
< A 11 c x + 2A 12 £ x«y + A 22 £ y + A 667xy 

V 

* 2 * * 

+ Dffw, xx + 2Df2 w >xx w >yy + ^16 w »xx w >xy (2) 

* 2 * * 2 

+ ^22 w »yy + ^26 w »yy w >xy + ^66 w »xy^ 

where c x , Cy, and y X y are the strain components. Commas denote differentiation with 
respect to the subscripted variables. 

The strains are expressed in terms of the displacements u, v, and w using 
the nonlinear strain-displacement relations. 

In the derivations that follow, the coordinate variables x and y are 
normalized with respect to their respective panel dimensions. The normalized coordi- 
nates (£, t )) are given by 
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The strain energy in the frame is 
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The potential of the external loads is 

p p 

G = -bN xx u(l,? 7 ) dr} - bN xy v(l,Ty) dr} 


( 6 ) 


The solution method employs the principle of minimum potential energy. In 
applying the principle of minimum potential energy, the displacement components are 
assumed to be functions of the independent variables £ and ty . The selected functions 
must satisfy the displacement boundary conditions and minimize the total potential 
energy. A generalized series expression for the displacement functions with unknown 
coefficients was selected for the present analysis and they are as follows: 


u - + a l a £ 

v “ B nm f 2 + b l a £ (7) 

w ™ ^nm^3 + b nmf4 + w o^5 » 

The coefficients A^, B^, C^, a^ , and b^ are unknown coefficients to be deter- 

mined by minimizing the total potential energy. The term w 0 in Equation 7 is the 
initial imperfection at the panel center. The function — f 5 ( ^ , ^ 7 ) is the initial 
imperfection function in terms of the lateral displacement and satisfies the dis- 
placement boundary conditions. 

The total potential energy is minimized with respect to the unknown 
coefficients. The minimization process yields a system of nonlinear algebraic equa- 
tions. Details of these algebraic equations are given in Reference 10. These equa- 
tions can be expressed in the following form: 


dn 

3A 


L J 



( 8 ) 


where the subscript L denotes the linear terms of the partial derivative of the total 
potential energy with respect to a particular unknown coefficient (A) , subscript N 
denotes the nonlinear terms and C represents the terms that are independent of the 
unknown coefficients. 


The number of nonlinear algebraic equations in the system given by Equation 
8 depend on the number of buckling modes used in the analysis. The number of equa- 
tions can be calculated from the relation 4NM+2 , where N is the number of buckling 
modes in the x-direction and M is the number of buckling modes in the y-direction. 

As N and M increase, the number of nonlinear terms on the right-hand side of Equation 
8 also significantly increases. A large number of nonlinear terms present numerical 
difficulty in solving Equation 8 . On the other hand, although the postbuckling 
behavior of a stiffened panel is mixed-mode behavior in general, the displacement 
response is dominated by a single buckling mode. Therefore, if the dominant buckling 
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mode is known, the postbuckling behavior can be accurately described using a single- 
mode analysis. 

In addition, for the case of no initial imperfection, i.e., w Q - o, the 
total number of integrals involved in Equation 8 is reduced to 108. The reduced 
system of nonlinear equations can be solved with very high accuracy by an iterative 
technique using the method of successive linearization. In this method, each of the 
unknown coefficient, A^, B^, C^, D^, a^ , and is assigned an initial value and 
substituted into the right-hand- side of Equation 8. Equation now becomes a system of 
linear algebraic equations and can be easily solved for the new values of the unknown 
coefficients. Using the new set of coefficients as initial values, another set of 
improved coefficients can be obtained by solving the linearized system. This proce- 
dure is continued until the solution converges within a desired limit. In the actual 
solution, only the initial values of the coefficients at the first load level need to 
be assigned. At higher load levels, the initial values are obtained by extrapolating 
the converged solutions of the preceding load levels to reduce the number of itera- 
tions. This numerical procedure for single mode analysis was coded in the PACL1 
computer program. 

WEBSTER . This analysis is a generalization of the elasticity method used in 
References 8 and 12. The method models the stiffener flange and the adjacent struc- 
ture (skin) as separate orthotropic plates that are bonded together with a bondline 
of zero thickness. Each plate may have different mechanical properties. 

The local model considers a flange terminating at the skin, as shown in 
Figure 2. The problem is formulated as a generalized plane deformation problem 
(i.e., the stress and strain components are independent of the z - coordinates ) . 

The stress-strain relations are given by 

e i * a ik a k “ 1,2, 3, 4, 5, 6 (9) 


and 


- 0 if i * k for i,k * 4,5,6 

where a ^ are the components of the effective compliance matrix. The equilibrium 
equations are satisfied automatically by defining the stress functions 
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The equations governing the stress function (F) are derived from the stress-strain 
relations (Equation 9) and the linear strain displacement relations. The governing 
equation becomes 
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where 


0 i 1 - a t1 - tilhl , 1-1,2, 4, 5, 6 (12) 

a 33 

The homogeneous solution of Equation 11 can be approximated by a series 
expansion in terms of complex variables, defined as Z - x + /iy. The stress functions 
are assumed to be functions of the complex variables, specifically F — F(Z) . The 
governing equation becomes 

[022 + 4 2 (2£i2 + 066) + 4 4 01l] F iv (Z) - 0 (13) 

Thus, n, depending on material properties, are the roots of the characteristic equa- 
tion 


022 + A* 2 ( 2 012 + 066) + 4 4 011 “ 0 ( 14 ) 

The displacements can then be expressed in terms of the stress functions as 


2 

u “ (Mk 011 + 012) F '( z k) 

, 2 

v = — (Mk 012 + 022) F '( z k) 

4k 

where, (k = 1,2, 3, 4) is the kth root of Equation 14. 

The series expressions for the stress functions are then stated as 
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(16) 


where, are constants to be determined. The stresses and displacements can then be 
expressed in terms of the constants and the complex variables Z^. 

The constants aj^^) are determined by satisfying prescribed boundary 
conditions. The boundary conditions include zero stress along the free edges, and 
stress and displacement continuity conditions along the interface. These conditions 
are satisfied by expressing the stress and displacement components in polar coordi- 
nates (r,8). It should be noted that Equations 9 through 16 are derived for the skin 
and the flange, such that each one actually represents two equations. 

The following relations use the subscript (or superscript) 1 and 2 to denote 
the flange and the skin, respectively. Surfaces AB and FA are assumed to be stress 
free. These boundary conditions are stated respectively as 


( 1 ) 

°8 

( 2 ) 
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The bond between the flange and the skin is assumed to be ideal along the interface 
(i.e., the stresses and displacements along this line are continuous). The stress 
continuity conditions at the interface are stated as 

d) _ (2) 


°d 

(r , 0) 

- <70 (r ,0) 

(19) 

(1) 

(r , 0) 

(2) 


T re 

- t r 0 (r , 0) 

(20) 

The displacement continuity conditions at the interface are 


u(l> 

(r.0) 

- u< 2 )(r,0) 

(21) 

v(l) 

(r,0) 

- v( 2 )(r,0) 

(22) 


A system of eight homogeneous equations in terms of the eight unknown 
constants (i=l,2 and k=l,2,3,4) result when the stress and displacement expres- 

sions are substituted into the boundary condition relations, Equations 17 through 22. 

Global analysis of the structure provides the forces that are applied to 
surfaces CD and EF of the local region shown in Figure 2. Surfaces FA, AB, and BC 
are stress free. The overall applied force system must satisfy the equilibrium 
condition of the local region ABCDEF. WEBSTER solves this problem using boundary 
collocation with collocation points distributed along edges CD and EF. 

PACSTER . PACSTER is an integrated analysis package that utilizes the 
displacement field obtained from the PACL1 analysis to determine the interfacial 
stresses at the skin/stiffener junctions by performing WEBSTER analysis. Therefore, 
the stress formulation in WEBSTER must be modified. The displacements in the local 
(WEBSTER) coordinates system can be expressed in terms of the stress function (F(Z)) 
as given by Equation 15. This modification essentially changed the original pre- 
scribed boundary force problem of WEBSTER into a prescribed boundary displacement 
problem. The displacement field obtained from PACL1 analysis is then used as bounda- 
ry displacements in the WEBSTER analysis. This procedure has been automated and 
resulted in the present computer code PACSTER, which is operational on a 286/386 
computer . 


A flowchart detailing the logic of the analysis methodology is shown in 
Figure 3. The program requires material properties, panel geometry, sectional prop- 
erties of the stringers and the frames and the load parameters as input. Following 
input execution of PACL1 begins and the program calculates the displacement field of 
the entire panel. The displacements are then screened to identify critical loca- 
tions. Because the objective of this analysis is to determine the applied load at 
which skin/stiffener separation occurs, the critical locations are selected based on 
the maximum displacements in the vicinity of the stringer or the frame. A total of 
four critical locations are selected for WEBSTER analysis. They are selected based 
on critical v and w displacements near the stringers and critical u and w displace- 
ments near the frames. After the critical locations are identified, the program 
proceeds to transfer all the required data into the local coordinate system for the 
WEBSTER analysis. The displacements associated with the critical locations are used 
as boundary conditions for the interfacial stress calculations. WEBSTER is executed 
four times for the four critical locations identified. 
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The interlaminar normal and shear stresses obtained from the local analysis 
are then used in the failure analysis. The failure criterion currently used in 
PACSTER is an average stress criterion in a quadratic form. The average stress 
criterion suggested in Reference 13 is applied here. The average shear (r ave ) anci 
normal (o ave ) stresses over a characteristic length (a Q ) are computed, and failure is 
defined as 
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+ 

ave 



k T o , 




(23) 


where, r Q and a 0 are shear and normal strength of the interface material, respective - 

iy. 


The procedure described above is repeated for each load increment until the 
prescribed maximum load is reached. 


Figure 4 illustrates the typical critical location selection procedure. A 
curved panel enclosed with stiffeners and frames is considered. The stiffener flange 
is of length a s and the frame flange is of length bf. In-plane compression and shear 
forces are applied to the curved panel. The resulting displacement pattern obtained 
from PACL1 analysis is shown in Figure 4. For the skin/stiffener separation analy- 
ses, the region of significance is highlighted in the figure. From the displacement 
field, the program has identified S cz and F cz as critical locations. The out-of- 
plane stress distributions along the skin/stiffener interface at location S cz are 
shown in Figure 5. This figure indicates high stress concentrations, both normal and 
shear, occur at the junction of the skin and the stiffener. These stress concentra- 
tions contribute to the skin/stiffener separation failure of the panel. 


NUMERICAL EXAMPLE 


A curved panel tested in Reference 10 is selected as an example to 
illustrate the analysis procedure. The composite panel is 24 inches long, 12.2 
inches wide with a radius of 45 inches. The panel skin is made of A370-5H/3501-6 
woven graphite/epoxy and AS/3501-6 unidirectional graphite/epoxy. The material 
properties, skin layup, and stiffener configurations are shown in Figure 6. The 
flange of the frame has a layup of [452/902/0/902/4521- The panel is under axial 
compression with total load P x . The interfacial strengths are assumed to be o o -3000 
psi and r o “11000 psi. These allowables are typical of the values obtained from 
flatwise tension and short beam shear tests, respectively. 

The results of the PACSTER analysis are shown in Figure 7 . The critical 
locations identified by the program are shown in the upper- left corner of the figure. 
They are denoted by (1) v-s/s: critical skin/stiffener interface due to v-dis- 

placement; (2) w-s/s: critical skin/stiffener interface due to w-displacement ; (3) 

u-s/f: critical skin/frame interface due to u- displacement; and (4) w-s/f critical 
skin/frame interface due to w-displacement. As shown in the figure, these locations 
depend on the applied load level. 

The interfacial stresses are expressed in terms of the failure index (Fj) 
defined in Equation 23. Figure 7 shows that the most critical location is w-s/s. At 
this location, Fj = 1.0 when the total compression reaches 9100 lb. That is 
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skin/stiffener separation failure will be initiated at the location w-s/s at 9100 lb 
of axial compression. Test results in Reference 10 showed the initial buckling 
occurred at 325 lb/in or approximately 4000 lb, and final failure of the panel, by 
skin/stiffener separation, at 825 lb/in or 10,050 lb. The results in Figure 7 show 
that PACSTER under-predicts the final failure by approximately 10 percent. This may 
be due to the conservative estimate of the interfacial strengths (a Q and r Q ) or the 
uncertainty in the out-of -plane properties of the materials. 


SUMMARY 


A methodology has been developed to predict the failure of postbuckled, 
stiffened composite panels by skin/stiffener separation. The methodology is applica- 
ble to curved as well as flat panels under combined uniaxial compression and in-plane 
shear loads. The entire analysis can be performed on a personal computer. 

Limited correlation of analytical predictions with experimental data 
indicated that the analysis provided reasonable prediction of failure load. Further 
data correlation is necessary to validate the methodology. 
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Figure 3. Flowchart for Analysis of Skin-Stiffener Separation in 
Postbuckled Panels. 



Figure 4. Postbuckled Panel With Stiffener Termination Details. 













Figure 5, Interlaminar Stress Distribution at Stiffener Termination 
in Stiffener Critical Zone. 


Panel Configuralion: Hat stiffened, 45-in radius, 24-in length, 12.2-in 

stiffener spacing 

Material A370-5H/3501 -6 woven graph rte/epoxy 

AS/3501-6 unidirectional graphite/epoxy 

Lamina Properties 



A370-5H/3501 -6 

AS/3501 -6 

Per ply thickness, in 

0 013 

0.0052 

EL. psi 

10.0 x 10® 

18 7* 10® 

EC PS‘ 

9.2 x 10® 

1 87 x 10® 

GLT. psi 

0.9 * 10® 

0 85 x 10® 

NULT 

0 055 

0.3 



Compression Uliimate Strains e"" - 0 012 lor layups with at least 

40% 0 degree plies 

Ec" - 0 015 


An underscore denotes woven ply 


Figure 6. Example Problem Data for Program PACSTER 











